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Abstract 

This is a sequel to a previous detailed study of quantum corrections to cos- 
mological correlations. It was found there that except in special cases these 
corrections depend on the whole history of inflation, not just on the behav- 
ior of fields at horizon exit. It is shown here that at least in perturbation 
theory these corrections can nevertheless not be proportional to positive 
powers of the Robertson- Walker scale factor, but only at most to powers of 
its logarithm, and are therefore never large. 
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I. INTRODUCTION 



Calculations of non-Gaussian corrections to cosmological correlations in 
the classical approximation have shown that these higher-order corrections 
are generally suppressed by powers of GH 2 , where G is Newton's constant 
and H is the cosmological expansion rate at the time of horizon exit[l]. 
From the magnitude of the observed Gaussian correlations of fluctuations 
in the cosmic microwave background, it is known that GH 2 ~ 10~ 12 for 
the fluctuation wavelengths studied in the microwave background, which 
makes the expected non-Gaussian corrections quite small. Quantum effects 
involve additional powers of G, and are therefore usually supposed to be 
too small to be detected. But a recent paper[2] has shown that it in many 
theories these quantum corrections depend on the whole history of inflation, 
not just on the behavior of fields at horizon exit, raising the possibility that 
they may be much larger than usually thought. In particular, if quantum 
corrections were to involve positive powers of the Robertson- Walker scale 
factor a(t) at the end of inflation, then they might be large enough to be 
detected. In this paper we will extend the results of reference [2] to a very 
large class of theories, and show that (at least in perturbation theory) this 
never happens; quantum corrections depend at most on powers of lna(i), 
and therefore (without sa 10 12 e-foldings after horizon exit) never become 
large. 

II. SCALARS AND GRAVITATION 

We will first consider a theory of multiple scalar fields (f n (x) and gravita- 
tion, slightly more general than that considered in reference [2]. The scalar 
Lagrangian is assumed to consist of a conventional minimal kinematic term, 
plus a term with arbitrary potential V(<p). In the Arnowitt-Deser-Misner 
formalism[3], the components of the metric are 

gij = a 2 e 2C [exp7]y, Ju = , (!) 

5 oo = -N 2 + g t3 N % Ni , g i0 = 9ij W , (2) 

where a(t) is the Robertson- Walker scale factor, 7y(x, t) is a gravitational 
wave amplitude, £(x, t) is a scalar, and N and N l are auxiliary fields, whose 
time-derivatives do not appear in the action. The Lagrangian density (with 
8ttG = 1) is 
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where Vj is the three-dimensional covariant derivative calculated with the 
three-metric e 2 ^.,-; and is the curvature tensor calculated with this 
three-metric. The auxiliary fields N and N l are to be found by requiring 
that the action is stationary in these variables. This gives the constraint 
equations: 
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For our present purposes, all we need to know about R 
solutions for N and N' 1 is that none of them contain terms with positive 
powers of a. We can impose gauge conditions, for instance by setting any 
one of the scalar fields equal to its unperturbed value and requiring that 

Qiiij = 0. 

The possibility of positive powers of a(t) in correlation functions of the ip n 
and/or £ arises from the explicit factors a 3 and a in the terms in Eq. (3). But 
these factors can be compensated by negative powers of a(t) in various field 
time derivatives and in various commutators, which arise from the structure 
of the perturbative expansion for correlation functions. The expectation 
value of any product Q{t) of field operators at various space points (but all 
at the same time t) is [2] 
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(with the N = term understood to be just Here Q 1 is the 

product Q in the interaction picture (with time-dependence generated by 
the part of the Hamiltonian that is quadratic in fluctuations); and Hj is 
the interaction part of the Hamiltonian (the part that is of third or higher 
order in fluctuations) in the interaction picture. The fields in the interaction 
picture are 



C(x,t) = J d 3 9 [e^ x a(q)C ? (t)+e-^ x a*(q)C,*(t) 



(8) 



^■(M) = f d 3 qJ2 [ eiq '%'(^ A)a(q, A) 7 ,(t) + e^e^g, A)a*(q, A) 7 *(i) 



p n (M) = J d 3 g[e iq - x a(q,n)^(t)+e- iq - x a*(q,n)^(t) 



(9) 
(10) 



where A = ±2 is a helicity index and eij(q, A) is a polarization tensor, while 
a(q), a(q, A), and a(q, n) are conventionally normalized annihilation oper- 
ators, satisfying the usual commutation relations 



a(q) , a*(q') 



'(Q-q') 



a(q) , a(q') 
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a(q, A) , a*(q', A')] = <5 AA ,<5 3 (q - q') , [a(q, A) , a(q', A')] = , (12) 



and 



a(q,n), a*(q',n') = 5 n „/(5 3 (q - q') , a(q, n) , a(q', n') =0, (13) 

The expectation value in Eq. (7) is assumed to be taken in a "Bunch- 
Davies" vacuum annihilated by these annihilation operators. Also, C q {t), 
7q(t), and (p q (t) are suitably normalized positive- frequency solutions of the 
wave equations 
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The functions Cg(*)> 7g(*)) an d approach time- independent limits 

Cg, 7g, and 99° at late times during inflation, when the perturbations are 
far outside the horizon, with the remainders Cg(*) ~~ Cg; 7g(*) ~~ 7g> an d 
(t) — all vanishing essentially (apart from slowly varying quantities like 
H and e) as a~ 2 (t). In consequence, Cj(i), 7<?(i)> and <£>q(i) all vanish at 
late times like a~ 2 (t). Also, as shown in reference [2], the commutator of 
any two interaction-picture fields at times t\, £2 during inflation but long 
after horizon exit goes essentially as a~ 3 (t), with t either t\ or ti or some 
weighted average of the times between t\ and £2- The same is true for the 
commutator of a field and a field time derivative, while the commutator of 
two field time derivatives goes as a~ 5 (t). 1 

This asymptotic properties of the commutators can be seen by noting 
that 

[C(x,t) , C(y,Ol = / A e^-yhm[c q (t)C q (t')} (17) 

[ 7 (x,i) , 7 (y,0] = / d 3 q e^( x -y)lm[ 7? (t)7*(i')] (18) 

b(x,t) , ¥»(y,0] = / A e^( x -y)lm[^(i)^(t')] • (19) 

The general solutions of Eqs. (14)-(16) are each linear combinations with 
complex coefficients of two independent real solutions, one of which goes at 
late times as a constant plus terms of order a -2 , while the other goes as a~ 3 , 
so the imaginary parts in Eqs. (17)— (19) arise only from the interference of the 
two independent real solutions, which goes as a -3 . Likewise the derivatives 
of these imaginary parts with respect to either t or t' also goes essentially 
as a~ 3 , because the derivative may act on the solution that already goes as 
a~ 3 , but the derivative with respect to both t and t' goes as a -5 , because 
both of the independent real solutions are differentiated. 

We are interested in the behavior of the correlation function (7) at late 
times t, when the perturbations are far outside the horizon. Inspection of the 

1 In this counting of powers of a(t), we are tacitly assuming that the time dependence 
can be evaluated before integrating over momenta, and will not be altered when the 
momentum integrals are done. This is based on the expectation that the counterterms 
introduced to eliminate ultraviolet divergences in flat space will suppress the contributions 
of large internal momenta even in an inflating spacetime. As discussed in reference [2], 
this expectation is not fulfilled for arbitrary choices of the operators whose correlation 
functions are to be calculated. It is necessary to consider only correlation functions of 
"renormalized" operators, for which large internal momenta em are suppressed. More 
work needs to be done to see how to construct appropriate renormalized operators. 
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Lagrangian density (3) shows that no term has more than 3 factors of a(t). 
According to Eq. (7), there are just as many commutators as interactions, 
and each commutator provides at least 3 factors of l/a{t) at late times, so 
the total number of factors of a(t) at late times in the integrals over time or 
in any subintegration is at most zero. With zero factors of a(t) the integrand 
can still grow like a power of t, which is more or less the same as a power 
of lna(t)[4], but it cannot grow like a power of a(t), and therefore (without 
~ 10 12 e-foldings) it cannot become large at late times. 

Indeed, since time derivatives of fields go like a _3 (t), and commutators 
of time derivatives of fields with each other go like a -5 (t) , the integrand will 
go like a negative power of a(t) if any interaction has less than 3 explicit 
factors of a(t), or if the time derivative of a field in any interaction does 
not appear in a commutator, or appears in a commutator with another time 
derivative. It is therefore only a very limited set of terms in the perturbation 
series that can contribute to a logarithmic growth of the integrand at late 
times. 

These conclusions would not be altered by the inclusion of higher-derivative 
terms in the action. Each pair of space derivatives is accompanied with a 
factor g l i oc a~ 2 , while field time derivatives of any order vanish at late 
times at least as fast as a~ 2 . 

It remains to consider the inclusion of other kinds of fields, but first we 
must say a word about the effect of scalar field masses. 

III. MASSES 

In the foregoing section we have treated the scalar fields (aside from a 
single inflaton field whose fluctuations can be eliminated by a gauge choice) 
as if they were all massless, with any possible scalar mass terms in the 
Lagrangian implicitly included as just additional possible terms in the po- 
tential V((p). As we have seen, when treated perturbatively such a term can 
at most introduce powers of In a in the late-time behavior of the integrand 
for cosmological correlation functions. But a mass m cannot be treated as 
a perturbation over time intervals t for which rat 3> 1, and in this case the 
powers of In a can add up to effects that materially change the late-time 
behavior of the integrand, requiring a separate treatment of mass effects. 

If a scalar mass m is sufficiently large compared with the expansion 
rate H, then it produces oscillations in the integrand at late times, which 
suppresses the contribution of any times later than 1/m. For m 3> H, the 
correlation function is therefore dominated by times in the era of horizon 
exit. But for m < H, a more detailed analysis is required. 
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We can get a good idea of what happens in these two cases by considering 
the simple example of a purely exponential expansion, a oc e Ht , with H 
constant. The wave equation for any one scalar field of mass m is 



m 2 + (q/a) 2 <p q = , 



For H constant, the solutions for q/aH <C 1 are 

2" 



C q a x ^ 



where 



l+0 (£)' 



+ D a a x ~ 




H 2 ' 



(20) 



(21) 



(22) 



and C q and D q are complex constants determined by matching this solution 
to solutions before horizon exit. For m > 3H/2 the exponents X± are 
complex conjugates, so as mentioned above, the oscillations of the wave 
functions suppress the contribution of late times. 
For m < 3H/2, the X± are real, with 

-3/2<A+<0, -3<A_<-3/2. 

Each scalar field factor in the Lagrangian thus contributes a factor of a A + (t) 
at late times, and as long as q/aH <C \J A+, the time derivative of a scalar 
field will contribute the same factor. On the other hand, commutators of 
scalar fields and/or scalar time derivatives contribute factors a(i) A++A ~ = 
a~ 3 (t), since the commutators can arise only from an interference between 
the two terms in Eq. (21). Once again, with no more than 3 powers of 
a(t) in each interaction, and with just as many commutators as there are 
interactions, the total number of factors of a(t) in the integrands for cor- 
relation functions cannot be greater than zero. Furthermore, except for 
trivial diagrams in which every vertex has just two lines attached, since 
each commutator involves just two fields, there must be fields that are not 
in commutators. These contribute additional factors of a(t) A+ to the inte- 
grand, and since A+ < 0, the integrand will be exponentially damped at late 
times, and the correlation functions will depend only on the behavior of the 
fields near horizon exit. 



IV. VECTOR FIELDS 
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Next consider a massless vector field, given (in temporal gauge) in the 
interaction picture by 

M*,t) = J2 [ A[e iq ' x £ 1 (g,A)a(q,AKW+e- ,t '' x 6 j (5,A) a *(q 1 A) H ;(t)] , 

A J 

(23) 

where here €i(q, A) and a(q, A) are the polarization vectors and annihilation 
operators for massless particles of helicity A = ±1, and u q (t) is a suitably 
normalized solution of the wave equation 

i(-(t)luS)) + ^u q (t) = (24) 

The commutator of two vector fields at unequal times is then 

[Ai(x, t) , A,-(x, t')\ = J d 3 q(5ij - qiq^j exp (iq ■ (x - x')) 

x(u g (t)u* q (t')-u q (t')u* q (t)) . (25) 

Now, the general solution of the wave equation (24) here takes the simple 
form 

u q {t) = C q cos qT + D q sin qr , (26) 
with C q and D q complex constants, and as usual 

f°° dt' , . 

TS l W)- (27) 

We see that at late times, where r — > 0, u q (t) approaches a constant, while 
u q (t) goes essentially as l/a(t). Also, 

u q {t) u*(t') - u q {t') u* q (t) = 2iIm(C 9J D*) sin (q(r - r')) (28) 

so at late times u q {t)u q (t') — u q (t')u q (t) and u q (t)u*(t') — Uq{t')u q {t) go es- 
sentially as 1/a, while u q {t)u* q {t') — u q (t')u*(t) goes to zero even faster, as 
1/a 3 . There are just as many commutators as there are interaction vertices, 
so if a term in the integrand involves a set of interactions H s with A s explicit 
factors of a, the integrand will contain altogether a number of factors of a 
bounded by 

#<E^" 1 ]- ( 29 ) 
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Because of the vector nature of the field, the maximum number of explicit 
factors of a(t) in any interaction is 3 — 2 = 1. For instance, in temporal 
gauge the electromagnetic interaction of a charged scalar field ip is 



a\~ 2 



ieAi(ip*di(p - ipdup*^ - e 2 AiAiip*(p , (30) 



with the factor a 3 coming from the metric determinant and the factor a~ 2 
coming from g %3 . So again the maximum number of factors of a in the 
integrand is zero, giving an integrand that grows at most like a power of 
In a. Derivative interactions of the vector field behave even better, because 
time-derivatives of vector fields give extra factors of 1/a, while pairs of space- 
derivatives are accompanied with factors oc a" 2 . For non-Abelian gauge 
fields A afl , there are self-interactions 



-a 3 a- 4 



C a p~jdiA a j AfoAjj + — C a f3yC a s e AipAj^A,isAj e 



(31) 



where C Q/ g 7 is a structure constant. The four factors of 1/a appear here 
because the interaction involves two contractions of space indices. Each 
such interaction contributes a factor a~ 2 to the integrand, suppressing the 
contribution of late times. 

V. DIRAC FIELDS 

A Dirac field of mass m in the interaction picture involves a wave function 
ip q (t) that satisfies the wave equation 

d 3H 

-J^q + + ^"SoTiftV'g + lomijjq = . (32) 

Hence for wave numbers far outside the horizon, the Dirac wave function 
has the asymptotic limit 

i/> q (t) cxe^ omt a- 3/2 {t) . (33) 

The matrix 70 has eigenvalues ±i, so the factor e~~ /omt produces an oscilla- 
tion, which does not affect bilinears like ipip or ?/'7oV ; > but does produce an 
oscillation in bilinears like "tpjitp, which suppresses the late-time contribution 
of interactions containing such bilinears. Even apart from this factor (as for 
instance for m = 0) , every bilinear combination of ip and ip is suppressed by 
a factor a~ 3 produced by the factor a -3 / 2 in Eq. (33). This in itself cancels 
the a 3 factor from the metric determinant, so that no positive powers of a(t) 
can be produced by any interaction involving Dirac fields. 
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VI. AFTERTHOUGHT 



In generic theories the N integrals over time in iV-th order perturbation 
theory will yield correlation functions at time t that grow as (hi a(t)) N . Such 
a power series in lna(t) can easily add up to a time dependence that grows 
like a power of a(t), or even more dramatically. As everyone knows, the 
series of powers of the logarithm of energy encountered in various flat-space 
theories such as quantum chromodynamics can be summed by the method of 
the renormalization group. It will be interesting to see if the power series in 
lna(t) encountered in calculating cosmological correlation functions at time 
t, though arising here in a very different way, can be summed by similar 
methods. 
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